We identify the asymmetry in reaction rates as a function of discharge/charge naturally within this approach.A molecular mechanical model for liquid water is developed that uses a physically motivated potential to represent Pauli repulsion and dispersion instead of the standard Lennard-Jones potential. The model has three atomic sites and a virtual site located on the ∠HOH bisector (i.e., a TIP4P-type model). Pauli-repulsive interactions are represented using a Buckingham-type exponential decay potential. Dispersion interactions are represented by both C6/r6 and C8/r8 terms. This higher order C8 dispersion term has been neglected by most force fields. The ForceBalance code was used to define parameters that optimally reproduce the experimental physical properties of liquid water. https://www.selleckchem.com/EGFR(HER).html The resulting model is in good agreement with the experimental density, dielectric constant, enthalpy of vaporization, isothermal compressibility, thermal expansion coefficient, diffusion coefficient, and radial distribution function. A graphical processing unit-accelerated implementation of this improved non-bonded potential can be employed in OpenMM without modification by using the CustomNonBondedForce feature. The efficient and automated parameterization of these non-bonded potentials provides a rational strategy to define a new molecular mechanical force field that treats repulsion and dispersion interactions more rigorously without major modifications to the existing simulation codes or a substantially larger computational cost.We analyze the probability distribution of rare first passage times corresponding to transitions between product and reactant states in a kinetic transition network. The mean first passage times and the corresponding rate constants are analyzed in detail for two model landscapes and the double funnel landscape corresponding to an atomic cluster. Evaluation schemes based on eigendecomposition and kinetic path sampling, which both allow access to the first passage time distribution, are benchmarked against mean first passage times calculated using graph transformation. Numerical precision issues severely limit the useful temperature range for eigendecomposition, but kinetic path sampling is capable of extending the first passage time analysis to lower temperatures, where the kinetics of interest constitute rare events. We then investigate the influence of free energy based state regrouping schemes for the underlying network. Alternative formulations of the effective transition rates for a given regrouping are compared in detail to determine their numerical stability and capability to reproduce the true kinetics, including recent coarse-graining approaches that preserve occupancy cross correlation functions. We find that appropriate regrouping of states under the simplest local equilibrium approximation can provide reduced transition networks with useful accuracy at somewhat lower temperatures. Finally, a method is provided to systematically interpolate between the local equilibrium approximation and exact intergroup dynamics. Spectral analysis is applied to each grouping of states, employing a moment-based mode selection criterion to produce a reduced state space, which does not require any spectral gap to exist, but reduces to gap-based coarse graining as a special case. Implementations of the developed methods are freely available online.Using a synchrotron-based Fourier-transform spectrometer, the high-resolution absorption spectra of the C1-symmetric 2,3-dihydrofuran (23DHF) and C2v-symmetric 2,5-dihydrofuran (25DHF) have been measured from 5.5 eV to 9.4 eV with an absolute absorption cross section scale. Oscillator strengths and vertical excitation energies of the lowest 18 states have been computed using the average of the second- and third-order algebraic diagrammatic construction polarization propagator method and the equation-of-motion coupled-cluster method at the level of singles and doubles model. These show that the bright valence transitions of ππ*-character are embedded into Rydberg transitions, whose oscillator strengths are at least one order of magnitude lower. To account for intensity borrowing, the first broad valence transition between 5.5 eV and 6.8 eV was simulated using a nuclear ensemble, and the agreement between experiment and theory is excellent. Whereas 23DHF only exhibits one broad valence transition followed by d/f Rydberg series converging to the ionization energy, the absorption spectrum of 25DHF has four bands, attributed to a valence nπσ → π*-transition, nπσ → 3px,z/3dxz transitions, a second valence nπ → π*-transition followed by d/f Rydberg series converging to the ionization energy, respectively. All Rydberg series converging to the ionization energy have been characterized in terms of their quantum defects.The excited-state properties and photoinduced charge-transfer (CT) kinetics in a series of symmetrical and asymmetrical Zn- and Au-ligated meso-meso-connected bacteriochlorin (BChl) complexes are studied computationally. BChl derivatives, which are excellent near-IR absorbing chromophores, are found to play a central role in bacterial photosynthetic reaction centers but are rarely used in artificial solar energy harvesting systems. The optical properties of chemically linked BChl complexes can be tuned by varying the linking group and involving different ligated metal ions. We investigate charge transfer in BChl dyads that are either directly linked or through a phenylene ring (1,4-phenylene) and which are ligating Zn or Au ions. The directly linked dyads with a nearly perpendicular arrangement of the BChl units bear markedly different properties than phenylene linked dyads. In addition, we find that the dielectric dependence of the intramolecular CT rate is very strong in neutral Zn-ligated dyads, whereas cationic Au-ligated dyads show negligible dielectric dependence of the CT rate. Rate constants of the photo induced CT process are calculated at the semiclassical Marcus level and are compared to fully quantum mechanical Fermi's golden rule based values. The rates are calculated using a screened range separated hybrid functional that offers a consistent framework for addressing environment polarization. We study solvated systems in two solvents of a low and a high scalar dielectric constant.